Constraints on supernovae dimming from photon-pseudo scalar coupling 
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An alternative mechanism that dims high redshift supernovae without cosmic acceleration utilizes 
an oscillation of photons into a pseudo-scalar particle during transit. Since angular diameter distance 
measures are immune to the loss of photons, this ambiguity in interpretation can be resolved by 
combining CMB acoustic peak measurements with the recent baryon oscillation detection in galaxy 
power spectra. This combination excludes a non-accelerating dark energy species at the 4(t level 
regardless of the level of the pseudo-scalar coupling. While solutions still exist with substantial 
non-cosmological dimming of supernovae, they may be tested with future improvement in baryon 
oscillation experiments. 

PACS numbers: draft 



I. INTRODUCTION 

Cosmic acceleration is mainly inferred through geo- 
metrical measures of the expansion history, most notably 
from the dimming of high redshift supernovae |ill,i3,SU|- 
In order to be certain of this interpretation, all other 
plausible explanations of the dimming must be ruled out. 

An alternative mechanism has been proposed for the 
dimming of supernovae 0. Here photons are converted 
into an unseen particle during transit. This mechanism is 
based on a generic interaction between the photon and a 
pseudo-scalar field. Other cosmological consequences of 
this pseudo-scalar (p-p) coupling include Faraday-like ro- 
tations, spontaneous generation of polarization, and pri- 
mordial magnetic fields generation |E S B 12 ■ 

Mediated by an external magnetic field, this coupling 
also produces the conversion that dims supernovae. This 
mechanism is similar to that exploited by axion direct 
detection experiments |lO| but in the reverse direction. 
Previous works examining the viability of this alternative 
dimming mechanism have largely focused on constraints 
from the observed achromaticity of the dimming |lllll2j |. 

In this Paper, we reconsider this alternate mechanism 
in light of recent cosmological measurements. First, we 
show that there exists large regions in the parameter 
space that both satisfy current luminosity distance mea- 
surements and avoid achromaticity constraints by satu- 
rating the conversion. Hence exclusion of this alternate 
explanation solely through supernovae data is difficult. 

On the other hand, cosmological probes of accelera- 
tion that do not involve luminosity distances are immune 
to dimming effects. Specifically, expansion history con- 
straints based on angular diameter distances from the 
measurement of standard rulers can be used to directly 
test the dimming hypothesis. The acoustic oscillations 
set by the photon-baryon fluid before recombination im- 



prints the sound horizon as a standard ruler in both the 
cosmic microwave background and galaxy power spectra. 
We combine measurements from WMAP which give 
the angular diameter distance to recombination with the 
recent detection of baryon oscillations in the SDSS LRG 
galaxy survey |14| which provides the distance ratio be- 
tween the galaxies and the CMB. 

Our result exclude a non-accelerating species of dark 
energy at the Aa level even allowing for the possibility of 
photon conversion. Furthermore, within the 95% confi- 
dence limits there exists separate solutions with no dim- 
ming and strong dimming of supernovae both of which 
require cosmic acceleration. These possibilities can be 
distinguished in the future with better constraints from 
the acoustic oscillations. 

In an Appendix, we also consider spectral constraints 
from COBE FIRAS on the dimming of the CMB it- 
self. These constraints, which are more robust than 
and equally powerful compared to CMB anisotropy con- 
straints , are useful for eliminating astrophysically dis- 
favored regions of the parameter space that involve large 
external magnetic fields with a short coherence length. 



II. PHOTON-PSEUDO SCALAR CONVERSION 

We study the interaction between photons and a 
pseudo-scalar field cj) through 
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where F^"" = \e^"'°'PFap and e'^""'^ is the Levi-Civita 
tensor. 

In a small patch where the effect of the cosmic expan- 
sion is negligible, the equation of motion for cj) is 
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W^cj) = -m^(j)- 
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-E-B. 
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The magnetic field B is a combination of the internal 
magnetic field of the radiation and any external mag- 
netic field that might exist. We assume however that 
within a domain the external magnetic field dominates. 
The external magnetic field in the intergalactic medium 
is expected to be around 10^^" Gauss with IMpc comov- 
ing domain size 0, 0| . 




FIG. 1: Top panel: the conversion probability with a given 
/i\/7 as a function of redshift z (see text). Middle panel: the 
magnitude change in the dimming between B and V bands. 
Bottom panel: dimmed bolometric magnitude of the super- 
novae in several cosmic acceleration scenarios parameterized 
by the dark energy equation of state w = p/p and the p-p cou- 
phng scenario. Data is plotted relative to the an accelerating 
ACDM model with w = —1 (horizontal line). 



We split the electric Maxwell equation into two com- 
ponents, parallel to B and E± perpendicular to B. 
The equation of motion for £'|| is 
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Here E\\ is coupled to the scalar field by the external 
magnetic field, but there is no mixture between E±^ and 
B. To summarize the coupled equations in matrix form. 
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where A — E/lo, ^ — B\\/M and I is the physical distance 
traveled in each domain. 

The off-diagonal interaction terms induce mixing be- 
tween A|| and 4>. The eigenvalues and the eigenstates of 



this mixing are given by 
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where the mixing angle is given by 



tan 26* 
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The conversion probability of a photon polarized in the 
direction in a single external magnetic domain is 



P-y_+^ = sin'^ 29 s\v? 



(8) 



where s is the physical size of a coherent field domain 
and ujp is the plasma frequency. As pointed out by [Tl| . 
the finite plasma frequency of the intergalactic medium 
sets a lower limit on the effective mass scale. We will 
hereafter assume that ujp ^ m throughout. 



III. SIMULATING DIMMING 

We now simulate the dimming of supernovae through 
multiple field domains in a cosmological context. We 
assume passive evolution of the parameters such that the 
field strength scales with the scale factor a as B (x a^, 
the domain size is fixed in comoving coordinates s cx a, 
and the plasma frequency 



1.2 X 10" 



fT^7^^)'''«"'/'eV, (9) 
V10-7cm-3/ ^ ' 



where Ueo is the free electron density today. 

A relevant parameter to consider is the ratio of the 
coherence length s to the plasma length Ip = 2u!/ cUp at 
the observation frequency [Tl|. If this ratio is large, then 
the spatially oscillating piece in the conversion proba- 
bility Eqn. will be averaged to 1/2 across the many 
random domains such that 
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in each domain. 



For optical light at t^sN 

-7r 



' IcV, the 

plasma length today Ipo ^ 100(neo/10~^cm~^)~^ kpc. 
With a typical magnetic domain size of order a Mpc, we 
expect s ^ Ip is a good assumption. Since the mixing an- 
gle scales with the plasma frequency as given in Eqn. l(7|l. 
we will take a fiducial value of = 2.4 x 10^^ cm~^ 
with the understanding that other cases may be obtained 
by replacing 



/i ^ ^ 
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Note that in this Hmit, the conversion probabiUty in 
a single domain is highly chromatic with a scaling of 

As an aside, there is also empirical evidence disfa- 
voring dimming in the opposite limit of s <C Ip. Here 
P-y^0 « 4s^/i^ and becomes independent of frequency. 
Hence low frequency photons from cosmological sources 
can be substantially dimmed. The blackbody nature of 
the CMB places significant constraints on dimming in 
this regime. In the Appendix we show that 



P-y~.<t>{uJcMB) < 2.5 X 10"^ (95%CL) 



(12) 



where wcmb = 1.24 x lO^'^eV. If at i^cmb the coher- 
ence length is still s <^ Ip then substantial dimming of 
supernovae would be completely ruled out. 

Under the assumption that the transition s ^ Ip lies 
somewhere between the optical and CMB frequencies, 
this constraint translates to one on Ipo / s through 



(30GHz vs 300GHz). The strong frequency scaling of 
P-y^0 cx Lu'^ more than compensates the lower preci- 
sion of the spectral measurements. With supernovae 
at zsN ~ 0.5 dimmed by order unity and assuming 
^CMB is high, e.g. comparable to the reionization redshift, 
s/lp > 1/10. Hence a very small coherence length to the 
field is not allowed if the supernovae are to be dimmed by 
order unity. We hereafter assume that s ^ Ip for viable 
models. 

To handle multiple domains, we employ a density ma- 
trix formulation. The initial conditions at emission are 
given by unpolarized photons, arbitrarily normalized to 
unity for convenience, and no excitation of the pseudo- 
scalar field 
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where D{z) is the comoving distance to redshift z, IpQ 
is evaluated at the SN frequency, and zcmb is the high- 
est redshift for which the magnetic field persists. Note 
that this 10~^ constraint from spectral distortions is 
more robust than and equally powerful compared with 
constraints based on the angular variation of the dim- 
ming. Although the latter lie at an amplitude of 10~^, 
anisotropy measurements are typically taken at a fre- 
quency that is up to an order of magnitude smaller 



In our notation for the density matrix 5f^(/), a and b 
denote the components (|| , _L, 0) with respect to the mag- 
netic field of the domain, the subscript i denotes the ith 
magnetic field domain, and the argument I denotes the 
physical location in the domain, for the beginning, Si 
for the end. 

Since the orientation of B in each domain is not iden- 
tical, the density matrix must be rotated between each 
domain. Defining ^ as the angle between B and the plane 
perpendicular to the propagation _B|| = cos^i?. With ip 
as the angle between in two neighboring domains, the 
density matrix transforms as 
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i sin 2il;{s}l^ - Si_\) + cos 2^S}^^ sin^ V'S'j'ii + cos^ i^S^k + sin 2i!S}\ sin ^I^S^ti + cos V'S'^^t | (^5) 
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where the right hand side is evaluated at position of the previous domain. We chose random variates uniform over 
27r for the angles in each domain. 

To propagate the density matrix across i, we employ the conversion probabilities implied by the mixing matrix 

Sf\s^i) = (1 - ^ sin^ 2e)sf\0) + i sin^ 293 + sin 20 cos 205,|''^(O) , 
S}'^{si) ^ i sin 29 cos 29 s} " (0) - ^ sin 261 cos 26*5/ "^(0) + sin^ 26'S'|'*(0) , 
Sf{s.i) = isin2205|"(O) + (l- isin220)5f'^(O)-sin20cos205|'^(O). 
There is no conversion in other components. 



The total conversion probability at the observer is de- 
termined by the appearance probability of the pseudo- 
scalar field at the last domain. The total conversion 



probability is given by 



■ 7^0 



i=D(z)f/S 
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where / is the one dimensional fraction of the external 
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FIG. 2: Qm — w constraints at the 68% and 95% CL. Cosmologies without p-p dimming are shown in the left panel and those 
with dimming in the right. Separate contours denote supernovae constraints (SN), baryon oscillations (R), and CMB acoustic 
peak constraints (CMB) as labeled. The ellipses at the intersection denote the combination of all constraints. 



magnetic domain between the observer and the source at 
a comoving distance D{z). Here S = Si/ai is the assumed 
constant coherence length in comoving coordinates. 

As can be seen in Eqn. the total conversion is 

a function of ii\/J . In the top panel of Fig ^ we show 
the conversion as a function of z for several choices of 
/iVT- Given a limit on the mass scale of the coupling of 
M > 2x lO^^eV [illlil, the magnetic field strength must 
be of order 10~^°Gauss to establish a substantial conver- 
sion across a cosmological distance. This is a plausible 
strength for the intergalactic field. 

The bottom panel of Fig ^ shows that the dimming 
produced by an accelerating component of dark energy 
with an equation of state w = p/ p < — 1/3 and by the p-p 
coupling with w > —1/3 looks almost identical (see the 
following section for model details). We would expect, 
as shown explicitly in the next section, that both scenar- 
ios would be equally favored by the luminosity distance 
measures alone. 

To mimic cosmic acceleration, the dimming must also 
be achromatic. Although the conversion probability in an 
individual domain is highly chromatic, the net dimming 
need not be so. The conversion probability equilibrates 
at 1/3 due to the oscillation of the pseudo-scalar back to 
a photon. Once the conversion probability is saturated, 
the dimming becomes achromatic. 

With a sufficiently large pVJ achromaticity can be 
achieved at a sufficiently low redshift to avoid obser- 
vational constraints. Observations limit the magni- 
tude difference of the dimming between B(0.44/iTO) and 
V(0.55/xm) bands to be at most 0.03 0. But as we 
see in the middle panel of Fig ^ once the conversion 



probability reaches the saturation limit, the difference in 
dimming between the bands disappears. 

Thus in terms of supernovae data alone, dimming by 
the p-p coupling and cosmic acceleration are equally fa- 
vored. 



IV. MODEL CONSTRAINTS 

We now explore the allowed parameter space of p-p 
dimming and cosmic acceleration models. The cosmic 
acceleration is parametrized by the dark energy equation 
of state w = p/ p and the non-relativistic energy density 
relative to critical flm in a flat universe. The p-p dim- 
ming model is parameterized by p-^/f. For simplicity, 
we will fix / — 2/3 and quote the allowed region in p 
evaluated at the present epoch. 

First, we use the gold set of the supernovae data alone 
0. When we consider the supernovae dimming just by 
the cosmic acceleration, a phantom dark energy model 
with w < — 1 is marginally favored as shown in the left 
panel of Fig|21 (see e.g. [13). It is interesting to note 
that the bolometric dimming of the supernovae most fa- 
vored by the cosmic acceleration scenario is statistically 
indistinguishable from the best fit p-p coupling scenario 
with w = —0.2 and no acceleration, as shown in the bot- 
tom panel of Fig ^ This degeneracy is quantified in the 
top-left panel in Fig|31 It shows that the cosmic accelera- 
tion by the phantom dark energy model with the negligi- 
ble p-p coupling is almost equally favored compared with 
strong p-p coupling with no cosmic acceleration i2(| . The 
two appear as discrete solutions with the ACDM model 
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in the weakly disfavored intermediate regime. 

Other cosmological tests are required to break this de- 
generacy. We next combine the luminosity distance mea- 
sured by the supernovae with the angular diameter dis- 
tance measured by CMB. The angular diameter distance 
at recombination z,, Da{z*) is given by fitting the CMB 
acoustic structures. We can estimate based upon 

the cosmological parameter ujm = ^mh^ = 0.135 ± 0.008 
which controls the physical extent of the sound horizon 
as well as the redshift of recombination and its measured 
angular scale = 301 ± 1 Ellp. 

The right panel in Fig [21 and the bottom left panel in 
Fig 121 show that the p-p coupling scenario with no cos- 
mic acceleration is now more likely than either phantom 
models or ACDM a.t w ~ —1 but the latter is still allowed 
at 95% confidence. 
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FIG. 3: w- fi constraints at the 68%, 95% and 99% CL. The 
magnetic field filling fraction / is fixed to 2/3 for simplicity. 



The favored p-p coupling with the current dataset 
requires a strong external magnetic field with ^^/J ~ 
8x10"'^° which saturates the conversion probability early 
in redshift. As shown in the middle panel in Fig^ the 
conversion probability with the strong external magnetic 
field passes the chromaticity constraint test. 



V. CONCLUSION 

We study an alternative explanation for the dimming 
of supernovae involving the loss of photons to a pseudo- 
scalar field (j> in transit. With only supernovae measure- 
ments on the bolometric total dimming and the achro- 
maticity of dimming, solutions without cosmic accelera- 
tion exist. 

This ambiguity in the supernovae data can be resolved 
with angular diameter distance measures that are based 
on a standard length scale instead of a standard candle. 
We apply the CMB acoustic peaks and the recently de- 
tected baryon oscillations as constraints on the model. 
We also use CMB spectral constraints to eliminate a 
class of small magnetic field coherence length solutions. 
We show that the non-accelerating solutions persist when 
adding just the CMB information but are eliminated at 
the 4(7 level once the baryon oscillations are added. 

Intriguingly solutions with substantial dimming ac- 
companied by cosmic acceleration still exist and could 
potentially bias measurements of the equation of state 
in the negative direction. These models can be tested 
with future baryon oscillation experiments. Likewise as- 
trophysical sources of non-cosmological dimming such as 
dust extinction can be tested by these means. 
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A much more direct test of p-p dimming comes from 
angular diameter distance measures to a redshift more 
comparable to the supernovae. This is provided by the 
recent detection of baryon oscillations in the SDSS LRG 
galaxy survey. This measurement can be summarized as 
a distance ratio 
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_ Da{z = 0.35) 
DAiz*) 



0.0979 ±0.0036. 



(17) 



When this constraint is added in, a non accelerating uni- 
verse with w > —1/3 is ruled out at the 4 a level (see 
Fig El). Interestingly, Fig. |3 shows two separate equally 
favored solutions: one with negligible dimming and the 
other with strong dimming. 



APPENDIX A: CMB SPECTRAL CONSTRAINT 

COBE FIRAS measured the CMB to be a pure black- 
body with no significant distortions A dimming of 
cosmological sources would distort the blackbody spec- 
trum and so in this Appendix we use the COBE FIRAS 
data shown in Fig. ^ to set limits on dimming scenarios. 

The p-p coupling dims the CMB spectrum as well as 
the supernovae light. We introduce the generic parame- 
ter Ai, representing the amplitude of the dimming of the 
CMB spectrum as 



1-^. ( 
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where the frequency v is measured in cm"^. Our lin- 
earized fit includes three unknown parameters, AT the 
temperature variation from To = 2.726K, Ay and the 
level of residual galaxy emission Gq . The form is 



1 



Vio 



where B^iTo) is the blackbody spectrum and §{1/) is a 
galaxy emission spectrum. We assume a fixed galaxy 
emission spectrum over the sky as given in |23l |. 

As shown in Fig.^ the best fit model is a pure black- 
body spectrum of Tg = 2.726K with a small high fre- 
quency galaxy emission contamination. However the 
shape of the linear distortion by a temperature uncer- 
tainty AT is mildly degenerate with A,^. Constraints on 
Ay are limited by this degeneracy as shown in Fig. ^ We 
obtain Ay = 1.0x 10-^±2.5 x f 0"^ at the 95% confidence 
level. 



FIG. 4; Residuals Io{i^) — Bi,{To) compared with model dis- 
tortions. Thin solid line denotes the best fit, thin dotted line 
denotes a deviation at the Ax^ = 1 level. This deviation 
is composed of compensating changes from the temperature 
variation AT, the p-p dimming A^, and the galactic compo- 
nent Go as shown. 
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